Stfeda-like formula in spin Hall effect 
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A generalized Stfeda formula is derived for the spin transport in spin-orbit coupled systems. 
As compared with the original Stfeda formula for charge transport, there is an extra contribution 
of the spin Hall conductance whenever the spin is not conserved. For recently studied systems 
with quantum spin Hall effect in which the ^-component spin is conserved, this extra contribution 
vanishes and the quantized value of spin Hall conductivity can be reproduced in the present approach. 
However, as spin is not conserved in general, this extra contribution can not be neglected, and the 
quantization is not exact. 
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Intrinsic spin Hall effect (SHE) offers new possibil- 
ity of designing semiconductor spintronic devices that 
do not require ferromagnetic elements or external mag- 
netic fields. This effect has been theoretically predicted 
both in p-doped semiconductors with Luttinger type of 
spin-orbit (SO) coupling and in n-doped semiconduc- 
tors with Rashba type of SO coupling pi* 2 ^ In the hole 
doped case^£ the transverse spin current is generated 
by the Berry curvature correction to the group velocity 
of a Bloch wave packet^* 5 - which is similar to the case 
of the charge current in quantum Hall effect £ Later it 
is pointed out that the intrinsic SHE can even exist in 
band insulators with SO coupling, which are called spin 
Hall insulators £ These spin Hall insulators would allow 
spin currents to be generated without dissipation. This 
dissipationless character of the intrinsic SHE again finds 
analogy in the quantum Hall effect. 

Pioneered by the work of Kane and Mele^ several 
specific single-particle Hamiltonians (graphena&SiiS and 
semiconductorsii*i2ii£) giving rise to the quantum SHE 
(QSHE) have been proposed, in which the intrinsic spin 
Hall conductance can be quantized in units of e/27r. 
These models can be considered as multiple copies of 
the charge Hall effect with different values of the spin, 
arranged so that the time-reversal symmetry is unbro- 
ken and the spin current is nonzero in the presence of 
an applied electric field. The investigations on its stabil- 
ity with respect to interactions and disorders have just 
begun, 10 : 14 : 15 

Because there exists similarities between SHE and 
quantum Hall effect, one may wonder if some concepts 
used in quantum Hall effect can be generalized to SHE. It 
is well known that (integer) quantum Hall effect can be 
analyzed by the Stfeda formula^ Stfeda showed that, 
if the Fermi level falls within the energy gap, the Hall 
conductance can be given by the charge-density response 
to a magnetic field (from orbital, rather than Zeeman 
coupling). This formula has been used to calculate the 
conductance of an electron gas in the presence of an ad- 
ditional periodic potential^ It has also been generalized 
to three dimensional systems^ Since the Stfeda formula 
is useful in quantum Hall effect, one may expect that 



similar formula for spin transport may be of some help 
in the study of SHE. Thus it is worthwhile to explore 
such a generalization^ 

In the present work, we show explicitly that, similar 
to the well-known result of Stfeda^ there are two terms 
in the (static) spin Hall conductance <j s h, one of which, 
a\ H , is due to the electron states at the Fermi energy, and 
the other one, cr^jy , is formally related to the contribution 
of all occupied electron states below the Fermi energy. 
Since the proposed models for QSHE thus far are all band 
insulators, where the density of states at the Fermi level is 
zero, therefore a^ H = 0. Hence we are mainly interested 
in the contribution to QSHE from crljj. Furthermore, we 
show that a^jj can be separated into a conserved part and 
a non-conserved part, in which the conserved part gives 
rise to a Stfeda-like contribution. That is, instead of 
directly calculating the spin-current response to an elec- 
tric field, we can calculate the spin-density response to a 
magnetic field (again from orbital, rather than Zeeman 
coupling) to obtain the conserved part of cr^jj. However, 
another contribution to a s ^j , which comes from the non- 
conservation of spin, is not zero in general. We note that 
our derivation is model independent, therefore it applies 
to any QSHE model. When this general formula is ap- 
plied to the aforementioned specific models of QSHE with 
conserved electron spin, the quantized values of o s h can 
be easily reproduced. As the spin in spin-orbit coupled 
system is not conserved in general, the non-Stfeda-like 
contribution in a J s jj can not be neglected, and it can re- 
sult in deviation from the quantized values, as shown in 
the numerical work of Ref. [lfl 

Using the linear response theory, the static spin Hall 
conductivity, for a z-component spin current flowing 
along the y direction under an electric field in the x di- 
rection, can be expressed as^°» 
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where il is the volume of the system, G ± (e) 
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lim^^ +( £ — H ± i-q)^ 1 is the operator Green function 
of electrons described by the Hamiltonian H. (j") 
is the charge (a-component spin) current density oper- 
ator in (i direction. The trace goes over every eigen- 
state in the space of the Hamiltonian H. In Eq. (JTJ, the 
delta function can be written in terms of Green functions, 
S(e -H) = - [G+(e) - G~(e)] /2m. 

If we keep half of Eq. and make an integration by 
parts on the second half, then a s n can be expressed as 
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Here a\ H stems from the contribution of electrons at the 
Fermi surface and a s jj formally contains the contribution 
of all filled states below the Fermi energy. Note that, at 



zero temperature, al H is proportional to the density of 
states at the Fermi energy. Therefore, for band insu- 
lators, where the Fermi level lies within the band gap 
such that the density of states at the Fermi level is zero, 
a\ H = 0. Since we are interested in the band insulators 
with QSHE, we need only consider the contribution from 

In order to simplify aj.jj , explicit expressions of j x and 
jy are needed. The charge current density operator is 
j x = —ev x /il. However, in studying SHE, the spin is in 
general not conserved due to the presence of SO interac- 
tion, thus the definition of the spin current is not unique. 
Here the conventional definition for spin current is used, 
that is, j* = {v y , s z }/2il. 

Using the cyclic property of the trace operation and 
the relations 
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we find that the spin Hall conductance a sH can be 
separated into a "conserved" part crllj^ 
conserved" part o^jy , 



and a "non- 
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def(e)Tv {s z G + (xv y - yv x )G + - s z G (xv y - yv x )G } , 
de/(e)Tr{s 2 (G+) 2 v x G+v y G + + G+v y G + v x (G+f - H. c. } 
j def(s)Tr{s z [G+[y,v x ] (G+) 2 + (G+) 2 [y,v x ]G+ +R. c] } . 
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Here we call a 



II, (n) 
H 



as the "non-conserved" part because it contains terms with s z = [s z ,H]/ih and therefore 



vanishes if s z is conserved. By employing the identity. 
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^{ Sz 6(s-H)} = £Tr{ Sz G+( 3 



the conserved part ct^h can be rewritten as a general- 

(9) 



ized Stfeda formula for SHE, 

II,(c) 
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where B is the magnitude of a uniform external magnetic 
field in the z direction, /i is the chemical potential, T is 
the temperature, and S z is the z-component spin density, 



(10) 



'->v ~ yv x )G + - s z G (xv y - yv x )G } , (8) 
I 

Eq. © shows the relation between the spin Hall conduc- 
tivity and the derivative of the z-component spin den- 
sity with respect to the perpendicular magnetic field B. 
Eqs. ©j 0, and © are the main results of this paper. 

Some comments are in order. Contrary to the case 
of the original Stfeda formula for charge transport^ 
there always exists an extra contribution for a^jj besides 
the Stfeda-like one as long as the spin is not conserved. 
Moreover, following the same reasoning of the present 
derivation, it is obvious that, for a current carrying some 
conserved quantity, there will always be a corresponding 
Stfeda-like formula. For example, when the total an- 
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gular momentum J z consisting of both spin and orbit 
angular momenta is conserved, the Hall conductivity for 
the total- angular-momentum current can be calculated 
through the total- angular-momentum density response to 
a magnetic field. That is, a Stfeda-like formula in which 
s z is replaced by J z can be derived. 

In the following, we illustrate the application of the 
present approach to some of the aforementioned QSHE 
models. As mentioned earlier, since the proposed models 
for QSHE thus far are all band insulators, thus what we 
need to consider is just the contribution from o~ s jj. It is 
shown in Ref. I 111 that, in the presence of a linear strain 
gradient, degenerate quantum Landau levels can be cre- 
ated by the SO coupling in conventional semiconductors, 
and then <j s h becomes quantized in units of e/2ir. This 
result can be understood as follows. Even though the 
intrinsic SO coupling is present, s z remains conserved in 
this case. Therefore, the non-conserved part c^^"' = 
and the spin Hall conductance comes from the Stfeda-like 
one [Eq. ©] only. To calculate the spin-density response 
to an external magnetic field B in the z direction, we 
notice that, at B = 0, opposite effective orbital magnetic 
fields ±B e g acting on spin f and J. electrons are induced 
through this special position-dependent SO interaction. 
When an external magnetic field B is turned on, the to- 
tal effective fields acting on spin f and J. electrons will 
become ±(B c s ± B). Therefore, for a fixed chemical po- 
tential /i, the number density of electrons with spin f and 
I becomes Ny± — e(B cS ± B)/h. Thus, using Eq. (jj), 

a S H = ^Yh — e/27r can be obtained. 

As the next application of our formula, we consider 
the QSHE model proposed in Ref. 0, where the intrin- 
sic SO coupling in a single-layer graphene film gives rise 
to QSHE with o~ s h being quantized in units of e/27r. 
This result can also be obtained by the present approach. 
For this model of graphene without the Rashba SO cou- 
pling, spin s z is conserved. Again, cr* H ^ = and the 
spin Hall conductance comes from the Stfeda-like one 
[Eq. ijjjjl] only. Because each independent subsystem of 
spin direction (f or J.) is equivalent to Haldane's spinless 



model, 22 following the same reasoning provided by Hal- 
dane, an extra field dependent number density of elec- 
trons for spin f and j will be AN^n = ±eB/h when 
an external magnetic field B is turned on. Thus, using 
Eq. (JSJ, a S H = = e /27r can again be reached. 

In the generic cases, spin is not conserved in the pres- 
ence of SO interaction and therefore spin current is not 
well-defined. As done in Ref. 0, one can separate the 
spin operator s z into a conserved part and a non- 
conserved part si™ : s z = si +s z n \ where the conserved 

(c) 

part s z consists of only intraband matrix elements of the 
spin. Thus a conserved spin current can be defined by 

(c) 

substituting the conserved part s z into the conventional 
expression of spin current. By employing this definition 
of spin current, it can be shown that the non-conserved 
part vanishes and the conserved part cr^ again 

obeys the Stfeda-like formula, Eq. , in which the spin 

operator s z is replaced by . In Ref. another defi- 
nition of the spin current is proposed which includes the 
torque dipole density. If we use this definition in the 
derivation for al H , the Stfeda-like formula Eq. © for 
the conserved part crf^y remains valid, and the non- 
conserved part cf^ n ' remains nonzero and has different 
expression from Eq. J7J. 

To conclude, a generalized Stfeda formula is derived 
for spin transport, and the applications of the present 
approach to the recently proposed QSHE systems are dis- 
cussed. Since the original Stfeda formula has succeeded 
in analyzing the quantum Hall effect in the presence a pe- 
riodic potential if the Fermi level lies within a miniband 
gap^ii^ we expect the generalized Stfeda formula will be 
of help for the spin transport under the same situation. 24 
Investigation along this direction is in progress. 
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